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Abstract. - We investigate the effect of quenched surface charge disorder on electrostatic 
interactions between two charged surfaces in the presence of dielectric inhomogeneities and 
added salt. We show that in the linear weak-coupling regime (i.e., by including mean-field 
and Gaussian-fluctuations contributions), the image-charge effects lead to a non-zero disorder- 
induced interaction free energy between two surfaces of equal mean charge that can be repulsive 
or attractive depending on the dielectric mismatch across the bounding surfaces and the exact 
location of the disordered charge distribution. 



Electrostatic interactions are one of the two fundamental components of the DLVO theory 
of colloidal stability [1,2]. They are standardly described by the Poisson-Boltzmann (PB) 
theory [1, 3] embodying the mean-field approach to classical charged systems. Mean- field 
interactions between like-charged macroions are repulsive in nature and thereby tend to sta- 
bilize solutions of charged macroions. In strongly coupled systems (e.g., when multivalent 
counterions are present), electrostatic interactions however induce strong attractive forces be- 
tween like-charged macroions [4,5], and thus act more like Lifshitz-van der Waals interactions 
that tend to destabilize charged solutions. This attraction can not be captured by the mean- 
field approach and a new paradigm dubbed the strong- coupling limit [6, 7] was introduced 
to describe the equilibrium properties of Coulomb fluids when the mobile counterion charges 
become large. The crossover from the mean-field Poisson-Boltzmann description to the strong- 
coupling limit is governed by a single dimensionless electrostatic coupling parameter, which 
is given by the ratio of the Bjerrum length (identifying Coulombic interaction between ions 
themselves) and the Gouy-Chapman length (describing electrostatic interaction between ions 
and the charged macroion surface) [4]. Electrostatic interactions between charged macroions 
in the mean-field and the strong-coupling limit thus unfold into a much richer structure than 
conveyed for many years by the DLVO paradigm. The collapse of a highly charged polyelec- 
trolyte, such as DNA, in the presence of multivalent counterions is the most dramatic example 
of unexpected and counter-intuitive features of the strong-coupling electrostatics [7,8]. 

Recently we added a new twist to the theory of electrostatic interactions in charged sys- 
tems [9]: not only can electrostatic interactions between like-charged macroions turn from 
repulsive to attractive due to strong-coupling counterion-induced correlations, but we showed 
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that in the case of no added salt and no image interactions, the quenched disordered distri- 
bution of surface charges on the macroions can induce an additive attractive interaction in 
the strong-coupling limit even if the mean charge of the macroions is zero. This effect is due 
to the nonlinear features of the average over quenched disorder of the distribution of charges 
on the macroion surfaces. Such quenched distributions of macroion charge have been invoked 
recently in experimental investigations of interactions between solid surfaces in the presence of 
charged surfactants [10, 11]. The patterning of interacting surfaces by quaternary ammonium 
surfactants in these experiments is highly disordered, depends on the method of preparation 
and has basic implications also for the forces that act between other types of hydrophilic 
surfaces with mixed charges, most notably in biological as well as in synthetic systems. Moti- 
vated by these observations we will now try to develop the theory of electrostatic interactions 
in systems with quenched disordered macroion charge distributions further. 

Here we will consider the effects of added salt and of image interactions due to dielectric 
inhomogeneities on the disorder-induced interaction between two charged walls of equal mean 
charge density. We shall focus only on the weak-coupling regime and evaluate interaction 
free energies up to the first-loop (Gaussian-fluctuations) contribution around the linearized 
mean-field (Debye-Hiickel) solution. We will show that in general image interactions have a 
pronounced effect on the way disordered charge distributions bring about electrostatic inter- 
actions in salt solution, their most notable effect being that they can induce non-monotonic 
interactions as a function of the spacing between the interacting surfaces. These results 
markedly contrast the weak-coupling results obtained in the absence of added salt and image- 
charge effects [9], where the disorder contribution turns out to be nil. 



Fig. 1 - Geometry of a system composed of two surfaces with disordered charge distribution at 
separation D = 2a with different dielectric constants for the interior region (where salt ions are 
present), e m , and for the exterior region, e p . 

Assume a monovalent salt solution confined between two charged surfaces at z = ±a. As- 
sume furthermore that the surface charges are small such that the linearization approximation 
is valid. The corresponding grand canonical partition function in the field of an external fixed 
charge distribution, p(r), is given in the functional integral representation as 




-m 




(i) 
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where </>(r) is a fluctuating potential field and the linearized effective Hamiltonian reads [12] 

S[<f>(r)] = ^ J dr e(r) (V0(r)) 2 + K 2 (r)0 2 (r) + J dr p(r)0(r). (2) 

Here we take the dielectric constant e = e(r) to have a value, e p , in the two semi-infinite regions 
(i.e. \z\ > a), which is in general different from, e m , assumed in between the charged surfaces. 
Accordingly, the inverse Debye screening length is n(r) = n = \Z8n£-Bn in between the two 
surfaces (with £b = e 2 /(47reoefcBT) being the Bjerrum length and n the salt concentration), 
and zero otherwise. We assume that the surface charge distribution, p(r), has a quenched 
disordered component. The average over quenched disorder is done as 



F=-k B T\nZ, (3) 
where the disorder average is defined here via a Gaussian probability distribution as 

(70= fv[p(r)](...)e-U dr ^ 1 (r)(p(r)- P o(r))\ (4) 



Note that po( r ) represents the mean charge density and g(r) gives the disorder variance around 
the mean value. Since all the functional integrals in the expression for the disorder-averaged 
free energy are Gaussian, the free energy follows straightforwardly as 

T = -k B T\nZ = | Tr 5 (r)g(r,r') + ^ Tr In Q~ x (r, r') + \ I Idvdr' p (r)£(r, r> (r')- (5) 

Here we have defined the inverse of the operator (3£eo(— V 2 + n 2 ) as the Green function that 
satisfies 

/3 £m£o (-V 2 + K 2 )g(r,r')-<5 3 (r-r') (6) 

with the appropriate boundary conditions of the continuity of derivatives multiplied by the di- 
electric constants at the surfaces with dielectric discontinuity. The disorder-averaged partition 
function could also be obtained through the replica formalism [9] but the direct integration 
approach is much more straightforward in the case of linearized effective Hamiltonian, Eq. (2) . 
In the second and third terms of Eq. w e recognize the usual fluctuational and linearized 
mean-field Debye-Hiickel (DH) contributions respectively. The first term is thus stemming 
from the effects of the disorder. Let us evaluate it explicitly and anlayze its consequences. 

Because of transverse isotropy, the following Fourier decomposition for the Green function 
is valid 

g(r,r') = J ^g(Q. z , z ') e -^<P-P'\ (7) 

with z and z' denoting the normal coordinate to the surfaces and p = (x,y), the transverse 
coordinates. We now evaluate two Green functions corresponding to the cases 

(i) when the disordered charge distribution is contained within the medium e m , and 

(ii) when the disordered charge distribution is contained within the medium e p . 

These Green functions can be derived straightforwardly by using the methods described in 
Ref. [13] as 



G(Q;z,z')= ' 



2/3e m e u 



2ae~ 4ua 

_ iua (e 2ua cosh u(z + z') + a cosh u(z - z')) 



(8) 
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for the case (i), and as 

G'(Q;z,z')= ' 



2/3e p e Q 
for the case (ii). Here 



— 4ua\ 

(9) 



-Q\z-z'\ _ "C 1 e ( e -Q\a+z'\+Q(a+z) 

1 - a 2 e~ iua V 



a{Q) ^ e m u{Q) e p Q ^ ^ = g2 + ^ 
e m M((y) + e p Q 

Furthermore, we assume that the disorder variance is surface distributed and thus 

g(r) = G S(z + a) + G 5(z-a). (11) 

It can reside either inside the slab (e.g., at a — S) or outside the slab (e.g., at a+5 for arbitrarily 
small 6). The only difference in the calculation is whether expressions (JHJ or @ are used when 
evaluating the first trace in Eq. [SJ Subsequently, one obtains the disorder-induced part of 
the free energy for the two aforementioned cases of the location of the disordered charge 
distribution as 

p T , w „ GS [°° n , n a (1 + a) 2 e- 4ua GSk„ . . .... 

for the case (i), and as 

Brr, , s~i, is GS [°° n a (1 - a 2 ) e- iua GS k . . 

for the case (ii), where S is the total area of the two bounding surfaces. In both of the 
above expressions we have subtracted the part of the free energy that does not depend on the 
separation a since we are only interested in the interaction free energy. 

It is thus immediately obvious that some asymmetry should exist in the system (either 
different dielectric constants in between and outside the surfaces, or salt in between and no 
salt outside, etc.) in order that the disorder contribution to the free energy becomes non-zero. 
Also obviously in both cases (i) and (ii) if a = 0, there is no disorder-induced interaction. If 
a = — 1, then in the case (i), the interaction is zero, but not in the case (ii). 

The dependence of the disorder part of the free energy, that is Eqs. ^] and ^3 on the 
dimensionless separation between the surfaces, na, is shown in Fig. |2 The most interesting 
feature of the disorder-induced interaction free energy is that the interaction can be non- 
monotonic and that it depends critically on the ratio of the two dielectric constants. For 
small and large values of na, the disorder interaction free energy assumes simple limiting 
forms as we show later. 

Putting the above results together with the fluctuational and mean-field contributions, 
where we assume that the surface charge distribution has a mean value given by 

po( r ) = c H z — a) + a 8(z + a), (14) 

we obtain the following expressions for the interaction free energy, Eq. [31 in the case (i) and 
(ii) respectively, i.e. 



T_ _ G f°° QdQ a{l + a) 2 e - iua , k B T 

S 47T£ m £ 



[°°QdQa 1 + a ) z e~ 4ua , k B T r° 2 _ 4ua . o* , , 

Jo u(l-a*e-^) jQWMl-c? e *»)+— (coth.a-1). 

(15) 
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Fig. 2 - Numerical evaluation of the disorder-induced interaction between two charged surfaces, 
F^(na) and i^hnffta) from Eqs. 1121 (left') and 1131 (right ) for < na < 1. The values of the ratio 
Em/£p are 0.2,0.4,0.6,0.8, 1.0,5.0 and 10.0 (from top to bottom). The non-monotonic character of 
the disorder-induced interaction is clearly discernible. Its details depend crucially on the ratio e m /e p . 



and 

P G [°°QdQ a (1 - a 2 ) e^ ua k B T 



r°° a 2 

/ QdQln(l-a 2 e~ 4ua )+ (cothwa-l). 

JO £m£oK 



S 47re p e J Q(l - a 2 e - iua ) 4tt 

(16) 

These are the final results of our calculation. We note here that the only approximation in- 
volved in the derivation of the above results is the linearization approximation in the Coulomb 
field action, Eq. [21 that makes them valid only in the weak-coupling limit, i.e., for small mean 
surface charge density a and low counterion valency. 

Let us assess the importance of the disorder-induced interaction by considering a few 
illuminating limiting cases, the general form being given numerically in Fig. [5] In the case of 
vanishing salt or small separations, na — > 0, one gets in the case (i), where disorder is located 
inside the slab of dielectric constant e m , 

{n rf, )g (.,0- g|^^ /(^). (17) 

where 

r°° du e' u 

'<■>-.£ Ti^r (18) 

The above limiting form is valid only if e m ^ e p . Otherwise higher order terms come into 
play. The second case (with the disorder located outside the slab and in the medium e p ) leads 
to exactly the same free energy and thus in this limit, there is no difference in the disorder- 
induced interaction whether the disordered charge distribution is within medium e m or e p . 
Obviously in this limit, the disorder- induced part of the interaction falls off inversely with the 
first power of the separation, D — 2a, to be compared with the inverse-square decay in the 
case of the zero-frequency van-der-Waals (fluctuational) term. Its sign depends on the values 
of both dielectric constants. One should also note that in this limit, the disorder and the 
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mean-field term combine, yielding 



■T 7 _ c 2 | G £ m (e m — e p ; 



S" £ m £oK 2 a \ 4:TT(T 2 (e m + £ P ) 3 \£ m +£ P JJ 

It would thus seem that the disorder merely renormalizes the square of the charge density. But 
since the disorder term can be either positive or negative, depending on the value of e m — e p , one 
can not claim that the only effect of disorder in this limit is the disorder-renormalization of the 
mean surface charge, since the whole expression Eq. 1191 can not be written as proportional 
to {<Jr) 2 = (<J — c') 2 , which is by definition always positive, where a' would indicate the 
disorder dependent terms. The disorder in this limit therefore does not simply renormalize 
the surface charge and can lead to attractive or repulsive interactions, depending on the sign 
of e m — e p . However, an important consequence of the disorder effects in this limit is that it 
induces interactions even between nominally uncharged surfaces with a mean charge density 
(7 = 0. These interactions have the same dependence on the separation as the mean-field DH 
term in this limit, except that they can be cither repulsive or attractive depending again on 
the difference e m — e p . Nominally neutral surfaces thus exhibit electrostatic- like interactions 
induced solely by the variance of the charge distribution, not its mean value! 

In the opposite limit of large salt or large separations, na — > oo, one remains with 

S irsoe m 4a 16n (feci) s m eoK 

for the case (i). The disorder- induced component (first term) has the same separation de- 
pendance as the standard screened zero-frequency van-der-Waals (vdW) interaction (second 
term) , but is shorter ranged than the corresponding mean-field DH term (third term) . Also in 
this limit, the disorder-induced interaction is always repulsive, which means that the overall 
interaction can change sign upon increase of the separation, as is already apparent from Fig. 
121 The interesting point now is that the disorder-induced term clearly renormalizes the fluc- 
tuational (van-der-Waals) contribution, since it has the same separation-dependance as the 
zero-frequency van-der-Waals term but with the opposite sign. 

For the second case (ii) and in the same limit of na — > oo, we obtain 

S ~K£Q£ m V & Ka 16-7T (na) e m £o K 

Again the disorder part of the interaction (first term) has almost the same functional de- 
pendence on the intersurface separation as the zero-frequency van-der-Waals part (second 
term). 

One can thus make a general conclusion that in the limit na — > 0, the disorder- induced 
component of the interaction free energy effectively behaves like the mean-field contribution, 
while in the limit na oo, it behaves like the fluctuational (vdW) contribution. In a certain 
sense, the disorder-induced interaction thus interpolates between mean-field and fluctuational 
interactions. All this is of course valid only in the weak-coupling limit and one can not apply 
these conclusions to the disordered strong-coupling regime [9] . The analysis of the interplay 
between disorder-induced effects and image-charge effects in the strong-coupling limit will be 
left for a separate exercise. 

The results derived above, apart from the effects due to finite salt concentration in be- 
tween the apposed charged surfaces, clearly differ from those obtained in the absence of 
dielectric inhomogeneities [9], where the mean-field contribution was shown to be decoupled 
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from the disorder effects. The conclusion reached in Ref. [9] is thus limited to disordered 
charge distributions immersed in a single dielectric medium without any inhomogeneities in 
the corresponding static dielectric constant. 

Note that if the mean surface charge is zero (a — 0), then according to Eqs. 1201 and 
1211 the equilibrium spacing between surfaces is given by the competing disorder-induced and 
fluctuational interaction parts for large na. In case (i), Eq. 1201 obviously implies no finite 
equilibrium spacing. The interaction is monotonic and its sign depends on whether the ratio 
4G/(K£ m £o&B? 1 ) is bigger or smaller than one. In case (ii), Eq. an optimal surface 
separation exists since the disorder and fluctuational contributions do not have exctly the 
same separation-dependence. 

In all the limiting cases addressed above, the disorder-induced part of the interaction can 
be masked by either mean-field or fluctuational terms in the total interaction, which would 
make its effects particularly difficult to pinpoint experimentally. Its most important feature, 
though, is the non-monotonic character of the interaction at intermediate separations (see 
Fig. [21 . In the case of interactions between charged interfaces, this feature may be important 
for the stability of planar charged macromolecular assemblies such as lipid bilayers. 
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